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ABSTRACT 

We report the detection of the Galactic nuclear disc in line-of-sight kinematics of stars, measured 
with infrared spectroscopy from APOGEE. This stellar component of the nuclear disc has an extent 
and rotation velocity V ~ 120kms _1 comparable to the gas disc in the central molecular zone. The 
current data suggest that this disc is kinematically cool and has a small vertical extent of order 50 pc. 
The stellar kinematics suggest a truncation radius/steep decline of the stellar disc at a galactocentric 
radius R ~ 150 pc, and provide tentative evidence for an overdensity at the position of the ring found 
in the molecular gas disc. 

Subject headings: stars: kinematics and dynamics - Galaxy: structure - bulge - center - disk - nucleus 


1. INTRODUCTION 

Through most of modern astronomy, empirical stud¬ 
ies of the central regions of the Milky Wav f ocused on 
radio observations of the central gas (cf. IBinnev et al.1 
119911 : iSormani et alll2015h . Visible photometry and spec¬ 
troscopy are nearly impossible due to the extinction from 
the Galactic mid-plane and the central regions them¬ 
selves. Eve n Baade’s Window i s at a latitude \b\ ^3 
degrees (cf. IGonzales et al.l [20125. While radio observa¬ 
tions have revealed the a ~ 200 pc gas disc in the central 
molecular zone (CMZ) of the Milky Way and its embed¬ 
ment in the Galactic bar both in density a nd gas kine¬ 
matics dRougoor fe Oortlll960t |Peterslll975D . stellar evi¬ 
dence is scarce. Apart from sporad ic stu d ies o n stellar 
kinematics via OH masers (see e.g. Habing et al.l F9831 
ILindavist et all 119921 IDeguchi et al.l 120041) . who found 
systematic rotation in the near-centre MASERS, the nu¬ 
clear disc has been indirectly inferred by modelling the 
photometric observations with estimated dens ity profiles 
(|Catchpole et al if 19901 ILaunhardt et al.l [20021) . 

The advance of infrared spectroscopy has now led to 
large systematic surveys of the Galactic mid-plane, in¬ 
cluding the Galactic centre, like the Apache Point Galac¬ 
tic Evolution Experiment (APOGEE, Majewski et al. in 
prep.). Earlier studies mostly focused on t he detectio n 
of sta rs around the central black hole (e.g. iSchodel et al.1 
l200l and using those orbits to co nstrain the black 
hole parameters an d distance (see iGhez et all 120091 : 
iGillessen et al.l[2009l) . Note that the nuclear disc studied 
here is a factor of 30 larger and likely of a different ori¬ 
gin than the ~ 5 pc di sc(s) around the central black hole 
(|Paumard et al.l [2006T) . Studies of the broader nuclear 
region of the Milky Way (see e.g. iSchultheis et, al.l [20031 
ICunha et al.llW0 71 were mostly limite d to rather small 
sampl e sizes or single clusters. E.g. IMatsunaga et al.l 
(120151) claim that the kinematics of their 4 Cepheid stars 
are consistent with the Nuclear Disc. 

The CMZ/nuclear gas disc is a direct consequence of 
the Galactic bar. Gas flows inward along the leading 
edge of the bar. Although some star formation has 
been detected along the leading edges of extragalactic 
bars dSheth et al.1 120021) . most of this gas will even¬ 



Fig. 1.— Overview of APOGEE stars (coloured dots) near the 
Galactic centre in Galactic longitude l and latitude b. Colours 
represent the mean line-of-sight velocity ui os of each star and its 
closest 29 neighbours. Note the division into plates/fields, and the 
clear dipole structure in ui os around the Galactic centre. 

tually c ross from a neighbou r ing x \-orbit onto the 
£ 2 -disc (Ivan Albada fe Sandersl 119821) . There it will 
accumulate until it reaches densities sufficient for star 
formation. Hence, some of the gas in the nuclear 
disc forms stars, a minor fraction is accreted onto the 
central black hole, and a major part will be expelled 
by stellar and Active Galactic Nucleus (AGN) feedback 
(cf. the bi polar outflow from the centre observed 
by iBland-Hawthorn fe Coheui 120031) . Although this 
central star formation has been studied in several 
pape rs (e.g. iSerabvn fc Morrisl 1 19961: Iv an Loon et al.1 
20031: iFiger et al.l 2004 lYusuf-Zadeh et all 120091 : 
Kruiissen et al.l 12014D . the nuclear stellar disc has not 
been directly detected from stellar spectra. Here we 
report the kinematic detection of the stellar nuclear disc 
using a sample of stars from the APOGEE survey. 

Data-related information is found in the next section. 
Section 3 presents observational evidence, a comparison 
to gas motions, and a qualitative exploration with toy 
models. We conclude in Section 4. 
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2. DATA 

We use stellar positions and line-of-sight velocities from 
APOGEE, downloaded from the Casjo bs system of the 
Sloan Digital Sky Survey III (SDSSIII, lEisenstein et al.1 
1201 If) . APOGEE is the first large infrared spectroscopic 
survey of the Milky Way. Its multi-object fibre spec¬ 
trograph simultaneously obtains about ~ 300 LL-band 
spectra with resolution R ~ 22500 per plate/exposure 
(see the plate structure in Fig. 1). The survey selection 
function (jZasowski et al.lf2012D aims at giant stars with 
the lowest possible complexity and bias. For each plate, 
all stars with (de-reddened) colour (J — Ks )o > 0.5 and 
de-reddened P/o-band magnitude within 6 < Hq < 11 
(for the central MW fiel ds) are selected with equal prob¬ 
ability. As discussed in IAume~ fc Schonrich l (j2015h . this 
selection favours young stellar populations (comparing 
ages of 1 and 10 Gyr) by about an order of magnitude, 
since the larger relative lifetime in evolved stages of mas¬ 
sive stars dominates over the opposite effect from the 
initial mass function. 

Wh ile the APOGEE selection function (Zasowsk i et all 
l2012f ) was designed to be unbiased, the Galactic Centre 
plate strongly favours positive Galactic longitudes. This 
asymmetry is not physical, but is caused by the survey for 
two reasons: manual selection and photometric crowding. 

The Ga lactic Centre field in A POGEE (termed ”GAL- 
CEN” in IZasowski et all 1201 2?) contains 233 out of 379 
stars that were added to the survey outside the normal 
selection criteria for comparison with pre-existing spec¬ 
troscopy. The good news is that neither APOGEE’S sec¬ 
ondary sample, nor the pre-existing samples were kine¬ 
matically selected. Hence, the detection of the disc kine¬ 
matics is fine, while the relative foreground contamina¬ 
tion and relative numbers within the structure may be 
distorted. 

The APOGEE ta rget s election uses 2MASS photome¬ 
try dSkrutskie et all 21)1)(») com bined with Rayleigh Jeans 
Colour Excess I Maiewski et al.1 1201 ll ) reddening esti¬ 
mates based on Spitzer-GLIM PSE data dBeniamin et al.1 
120031 : IChurchwell et all 120091) . The 2MASS photometry 
in this field is compromised by inhomogeneous extinction 
and crowding: the smaller extinction at negative longi¬ 
tudes, whic h is c on nected to the asymmetry of the CMZ 
(CMZ. iBallv et all 11988 1. and hence stronger crowding 
reduces the number of target candidates passing the qual¬ 
ity criteria. This prevents a clean assessment of the de¬ 
tailed structure of the nuclear disc. 

APOGEE does not offer stellar parameters for most 
objects in the Galactic centre. To diminish foreground 
contamination, we thus have to rely on photometric 
methods. Obviously, the handful of telluric standards 
can be removed by imposing (J—K )o > 0.5. We obtain a 
cleaner subsample via the lower extinction of foreground 
stars, requiring that stars have a valid reddening estimate 
> 3 based on WISE photometry fcf. iBlum et all 
119901 iFierer et, all 12004 [Gonzales et al.ll2012l ). This may 
preferentially remove stars on the front edge of the nu¬ 
clear disc, a bias which we assess by comparison of the 
two subsamples. 

To compare the stellar kinematics with the gas in 
the CMZ, we use publicly available data cubes from 
the AST/RO survey (Ant arctic Su bmillimeter Telescope 
and Remote Observatory, iMartin et al.112004 1). For the 


longitude-velocity (Iv) plot, we sum the emission temper¬ 
atures for all 73 latitude bins within the survey’s latitude 
range of —0.3 < b/ deg < 0.2 from the CO J = 4 —> 3 
lines at 461.041GHz. While this survey has a mildly 
low er resol ution compared to later data compilations (see 
e.g. Uones et al.l[20T2[ ). it offers the advantage of publicly 
available and well documented data, while it covers the 
main features observed in the Iv -plane by other surveys. 

We translate line-of-sight velocities ui os into Galactic 
rest frame velocities u ga i, defined as 

= v ios + U Q cos(Z) cos(6) + V Q sin(Z) cos(6) + W Q sin(6). 

( 1 ) 


For the solar galactocentric distance we assume Rq = 
8.27 kpc and for its motion against th e Galactic rest 


frame (LG, IG, W 

k) = (14, 250, 7) kms 1 (IGillessen et al. 

2009 

: McMillan 20111 Schonrich 2012t Schonrich et al. 

2010 

)■ 


3. KINEMATIC DETECTION OF THE NUCLEAR DISC 

The rotation of the nuclear stellar disc of the Milky 
Way is evident from Fig. [TJ limited to the innermost de¬ 
gree around the Galactic centre. Stars at negative lon¬ 
gitudes approach with mean ci os ~ —80kms _1 and re¬ 
cede with the same speed at positive longitudes. This 
pattern is limited to latitudes |6| < 0.4deg and is hence 
dir ectly connected to the di sky overdensity identified e.g. 
by lLaunhardt et alJ (|2002D . who find a disc radius of 
~ 230 pc and a scale height of ~ 45 pc (though the bulk 
of their 60/nn emission is within 1 degree of the Galactic 
centre). 

The case that this is indeed a rotating nuclear disc is 
simple: Despite the foreground contamination, the av¬ 
erage r>ios shows a shift of more than 100 km s -1 on this 
small scale. The feature is far too dominant for any back¬ 
ground effect behind the Galactic centre. Spiral structure 
may shape velocities by ~ lOkms -1 , an order of mag¬ 
nitude less than this observation. Also, its length scale 
would by far exceed the < 100 pc implied by the angu¬ 
lar size. The other candidate would be a con nection to 
the velocity features discovered by iNidever et al.1 (I2012D . 
However, these are again too weak in numbers and too 
extended covering an angular scale of A l > 10 degrees. 
Explaining a mean velocity shift centered around a mean 
< v\ os >~ 0 would require that miraculously the fore¬ 
ground feature is observed exclusively on one side and 
the background feature on the other side of the Galac¬ 
tic Centre. Last, no systematics are detected in the 
motion above |6| ~ 0.4 deg, i.e. altitudes larger than 
~ 50 pc, confirming a disc roughly within the Galactic 
plane. Hence, the only possible explanation is a rotating 
nuclear disc around the Galactic Centre. 

In Fig. [2] we show the velocity distributions for stars 
with |Z|,|6| < ldeg, dividing the field into stars with 
l > 0.05 deg and l < —0.05 deg. Already the entire sam¬ 
ple is dominated by the strong odd part in the veloc¬ 
ity distribution coming from the two sides of the disc 
at velocities around 50 — 120 km s. The small ve¬ 
locities are dominated by foreground stars in the disc 
mid-plane, while both sides show a couple of high ve¬ 
locity stars, which are associated with the bar/bulge 
region (stellar populations streaming alongside the bar 
have tails up to |twil ~ 300 km s -1 , see INidever et al.1 
l2012t lAumer &: Schonrichl 120151) . The disc kinematics 
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0.05 deg (208 stars) and l < —0.05 deg (129 stars) to present the two 
halves of the central disc. Even without any further selection, the 
two sides of the nuclear disc dominate the sample. The foreground 
contamination is reduced when we narrow the selection to highly 
reddened stars (A k > 3, 109 and 58 objects). 

are cleaned up via the extinction cut Ak > 3 (thicker 
lines) . 

We now compare the stellar kinematics in the Zu-plane 
to the motions of the molecular gas in the known nuclear 
gas disc. In Fig. [3] we display the integrated emission in 
the CO J = 4—1-3 transition with \b\ < 0.3 deg in u ga i 
(y-axis) versus longitude l with blue to white shades. We 
can clearly see the ridge of emission around 115 km s^ 1 
at positive longitudes, as well as a lower, sloping ridge 
connecting the bright emission of the Sgr B2 region at 
(t^gai) ~ (0.7deg, 70kms _1 ) to the origin. The lat¬ 
ter has been interpreted as a ring-shaped region of en¬ 
hanced radio emissio ns at a radius of R ~ 150 pc around 
the Galactic centre dMolinari et al.1 120111) . We overplot 
APOGEE stars with |f|, |&| < 1 deg (orange crosses). To 
curb foreground contamination again, we mark the high- 
reddening stars with Ak > 3 as red bullets. While the 
velocity distribution of low-reddening stars is nearly in¬ 
dependent of longitude, consistent with their being fore¬ 
ground stars, the high-reddening stars display the point 
symmetry typical for disc kinematics, plus moderate con¬ 
tamination. They also show some tentative structure 
along the two ridges defined by the molecular gas. 

Detailed quantitative exploration with consistent kine¬ 
matic models has to await a decent 3D-reddening map 
and reasonable stellar parameters to assess the fore¬ 
ground contamination. However, for a qualitative dis¬ 
cussion we provide three mock models in Fig. [Tj The 
plotted mock models are axisymmetric with a constant 
azimuthal streaming velocity V,$ ~ 120kms _1 . We shift 
the longitude of the toy model centre by A l = —0.054 deg 
to the position of Sgr A*. The data favour this, as the 
sign change in mean v ga \ occurs at negative longitude. 

The predicted velocity distributions are folded with a 
Gaussian kernel with a = 15kms“ to mimic the ran¬ 
dom motion of stars (which in reality should be neither 
isotropic, nor constant in R). The dispersion has an up¬ 
per limit by the sharpness of the sloping density ridge, 
and a lower limit by the expectation for secular heating. 
Model complexity increases from top to bottom, starting 
from a simple exponential disc. The middle panel intro- 


TABLE 1 

Fit parameters from equation(2) for stars with 5 < 0.4 deg. 
Errors are given by Aln(P) = l/v2 with full variation of the 
other parameters. 


name 

parameter 

all 

Ax > 3.0 

azimuthal vel. 

y<j>/ km s —1 

123.0±1, 9 

120lg 4 

scale length 

Rd/ pc 

52^17 

60 1™ 

radius of ring 

Rr/ pc 

150±lg 

147116 

ring/centr. density 

fr 

i-6±l;l 

3.0111 

trunc. radius 

Rt/ pc 

1501“ 

1471“ 

contam. fraction 

fc 

n 1 qc+0.055 
u.ioo_ 0 043 

0.059!°;^ 

contam. vel. disp. 

cr c / km s —1 

601J? 

54li3 

log likelihood 

HP) 

-440.2 

-325.4 


duces a truncation radius, and the bottom model adds 
a ring of increased density near the truncation. The toy 
model for the stellar surface density reads: 


E d (R) = e~ RIRd 


J 1 if R<Rt 

\ e -(R-R t )/R d ,2 ifR > Ru 


( 2 ) 


with cylindrical radius R, scale-length Rd, truncation ra¬ 
dius R t . We set the drop-off steepness Rd ,2 = lpc. To 
the bottom panel we add a 5 pc wide, constant surface 
density ring, centred on a radius R r and with density f r 
relative to the disc centre. The model without truncation 
has Rd, 2 —> 00 . 

The mild density maximum near the circular veloc¬ 
ity derives from the comparatively large distance range 
within which the azimuthal direction of the disc is nearly 
collinear with the line of sight. This is enhanced by the 
density profile, as the tangential parts probe the smallest 
radii. A radial truncation removes at fixed l the stars at 
small velocities, and limits the extent of the disc in lon¬ 
gitude. Note that the longitude limits of the APOGEE 
held do not affect our fits since we separately normalise 
narrow bins in longitude. A ring with enhanced density 
(see bottom panel) produces a straight density ridge be¬ 
tween the origin and the tangential longitude of the ring, 
where it ends in a pronounced density maximum. 

Fit parameters for the full mock model are listed in 
TableQ] In these fits we demand that R t < R r and al¬ 
low for a naive Gaussian foreground contamination with 
dispersion a c and a constant (in l) contamination frac¬ 
tion f c . Fits are restricted to positive l, sampled into 
0.1 degree wide longitude bins with separate normal¬ 
ization in star count. Negative longitudes are excluded 
due to the more uncertain selection function and any¬ 
way low number counts. We limit the fitting range to 
I'Cgai | < 130kms _1 . This curbs contamination from high 
velocity structures, i.e. f at ta ils in th e cont aminating ve¬ 
locity distributions (cf. lAumer fc Schonrichll2015l ). while 
still covering some of the decline in disc density above 
the rotation speed. Comparison between the results for 
samples with and without the reddening cut Ak > 3.0 
confirms that this cut removes most of the foreground 
contamination, while the fit parameters stay mostly un¬ 
affected. The data suggest a truncation radius, though 
its removal can, within the significance limits, be com¬ 
pensated by a shorter scale length (< 50 pc). If we drop 
the condition R t < R r the truncation radius moves in¬ 
wards to 80 pc with a significant improvement of the like¬ 
lihood. The fits demand a ring with comparable densities 
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Fig. 3.— Latitude-velocity ( Iv ) plot of stars and gas in the Galactic centre region. The colour code gives the intensity of CO J = 4 —> 3 
emission from the AST/RO survey integrated over all |6| < 0.3 deg in units of the maximum emission density. The ui os range not covered 
by the survey is left blank, the cut-off has a slope due to the translation into v ga \. Orange crosses indicate APOGEE stars at latitudes 
|6| < ldeg, red bullets the likely disc members with estimated A-band extinction A k > 3. 


to the disc centre (uncertain due to unknown width and 
dispersion), and with a radius around 150pc. Despite 
this seemingly robust result, we caution that these mock 
models are very unsatisfactory and should no be consid¬ 
ered more than a first guidance. The truncation radius 
is in a zone of heavy differential effects in reddening, and 
the APOGEE selection is not well constrained. 

These toy models underline the need for additional 
data: In particular the bottom panel of Fig. []] shows the 
need for a larger Galactic centre sample covering the sus¬ 
pected tips of the nuclear disc in longitude, which would 
also be important to assess a possible non-circularity. 


4. CONCLUSIONS 

We report the detection of the Galactic nuclear disc 
from stellar spectra of the APOGEE survey. 

Even without removing any foreground contamination, 
a dipole in line-of-sight velocities is prominent in the data 
with a velocity shift > 100 km s -1 , limited to within ~ 1 
degree from the Galactic Centre. The only viable expla¬ 
nation for such a strong mean velocity contrast on small¬ 
est scales is a rotating disc in the central ~ (150—180) pc. 

We show that the data are consistent with a nearly 
axisymmetric distribution with an azimuthal speed of ~ 
120 km s -1 . This is in line with gas motions of the CMZ 
derived from radio observations. Signatures of rotation i n 
the OH/IR objects were found bv ILindavist et all (1199211 . 
The disc appears confined to latitudes |6| < 0.4deg, i.e. 
an altitude \z\ < 50 pc c onsistent with the phot omet- 
ric density estimates from lLauuhardt et al.l (l2002lf . Due 
to the selection uncertainties and spatial limits, we do 
not use the spatial densities. Consequently the limited 
plate extent with |2| < 0.7deg does not alter our conclu¬ 
sions, but only limits our stated significance. The stellar 
kinematics alone provide evidence for a ring of enhanced 
density around R ~ 150 pc, outside which the data hint 
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Fig. 4.— Latitude-velocity plot of high-reddening APOGEE stars 
from Fig. [3] at |6| < 0.4 deg (red circles) compared with three toy 
models (blue shading) with increasing complexity: A simple ex¬ 
ponential disc (top), an exponential disc with a truncation radius 
beyond which the density drops steeply (middle), and a denser ring 
near the truncation radius (bottom, parameters from Table 1 for 
A k >3). In the top panel we add as green lines the line-of-sight 
velocities that perfect rings with rotational speed V = 120 km s~ 1 
and radii 50, 100, and 150 pc would have. Note that there are 
almost no APOGEE stars beyond |/| > 0.7deg due to the survey 
geometry. 


to a truncation/steeper decline. This ring feature co¬ 
incides with the molecular ring found in the gas mo¬ 
tions. Its small width in line-of-sight velocities suggests 
a kinematically cool population. The bulk of the nuclear 
disc mass resides within ~ 150 pc, with indications for a 
steeper decline outside the ring. Not allowing for a trun¬ 
cation, the scale length would drop to Ed < 50 pc with 
nearly the same effect. This implies a somewhat smaller 
radial scale than the ~ 230 pc found in photometry by 
lLaunhardt et al.l ( 2002 ), though it is in line with the bulk 
of their 60/mi emission. 

While the detection of the nuclear disc is beyond the 
tightest significance limits (without any clean-up the 
line-of-sight velocity signal from Fig. [T] exceeds 10cr), its 
precise nature is more ambiguous. We explore the pos¬ 
sible shape with toy models. However, any quantita¬ 
tive approach has to await the creation of 3D-reddening 
maps, and a manual analysis of the stellar parameters, 
to assess the precise foreground contamination and the 
selection function of stars in this region. On the other 
hand, this finding facilitates analysis of nuclear disc stars, 
as they can be clearly identified by their kinematics. 

These uncertainties prevent an assessment of possi¬ 
ble non-circularity in the nuclear disc. A mild non¬ 
circularity is expected, as this disc is born from gas on 
or bits of the x ? family and affected by the rotating bar 
ijBinnev et al.l I1991D . some might derive f rom possible 
spiral arms within this disc. IMolinari et al.l (120111) sug¬ 
gest an even more complex shape, though the kinematic 
stability of suc h a structure is doubtf ul; improvements 
are discussed in lKruiissen et alJ (120151) . _ 

In a recent paper. iDebattista et all (|2015l) claimed to 
have identified a ~ 1 k pc nuclear disc based on the high- 
velocity feature of iNidever et all ( 2012 ) at longitudes 
l ~ 10 degrees. Although the idea is intriguing, they 
have not presented conclusive evidence that their model 
can fit the obser ved r<i os distributions. The simulations 
of iSormani et all (|2015l) suggest such a large xi disc is 
unlikely. Also, it would require significant flattening of 
the density profi le sugg ested both by this study and by 
lLaunhardt et alJ (|2002f) . More critically, these obser¬ 
vations provide no p roof/detecti on of a ~ 1 kpc sized 
X ‘2 disc, since lAumer fe Schonrichl (120151 ) have identified 
them with orbits alongside the bar (largely x\ orbits), 
in a model that provides a good fit to the obser v ations . 
A similar conclusion was drawn by IMollov et all (I2015D . 
We would like to stress again that those high velocity 
kinematic features cannot explain the observations in the 
central degree: their stars are expected to have larger 
velocities (~ 200kms _1 ) than the nuclear disc, which is 
~ 100kms' 1 slower. More importantly they would be 
found on both sides of the Galactic centre, i.e. both the 
background stars at ~ +200 kms -1 and the foreground 
stars moving at negative velocity would be found on both 
positive and negative l. 

So far, APOGEE has no data available at longitudes 
at the expected tangent positions of the truncation ra¬ 
dius and the possible ring. Clearly, infrared surveys like 
APOGEE can penetrate into the nuclear disc and gather 
excellent data. A satisfactory study will require a 3D 
extinction map and samples covering the full range of 
longitudes at |Z| < 2 with a consistent selection function. 
Thus, this work should be understood as a starting point 
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for more comprehensive surveys of this exciting region. 
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